Spin-flavor oscillations of Dirac neutrinos described by relativistic quantum mechanics 
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Spin-flavor oscillations of Dirac neutrinos in matter and a magnetic field are studied using the 
method of relativistic quantum mechanics. Using the exact solution of the wave equation for a 
massive neutrino, taking into account external fields, the effective Hamiltonian governing neutrino 
spin-flavor oscillations is derived. Then the The consistency of our approach with the commonly 
used quantum mechanical method is demonstrated. The obtained correction to the usual effective 
Hamiltonian results in the appearance of the new resonance in neutrino oscillations. Applications to 
spin-flavor neutrino oscillations in an expanding envelope of a supernova are discussed. In particular, 
transitions between right-polarized electron neutrinos and additional sterile neutrinos are studied 
for realistic background matter and magnetic field distributions. The influence of other factors such 
as the longitudinal magnetic field, the matter polarization, and the non-standard contributions to 
the neutrino effective potential, is also analyzed. 
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I. INTRODUCTION 

It was confirmed by numerous experiments that neu- 
trinos are massive particles and there is a mixing be- 
tween different neutrino generations. Besides these non- 
standard model neutrino properties it is believed that 
neutrinos can also have non-zero diagonal and transi- 
tion magnetic moments. The latter mix both eigenstates 
with opposite helicities as well as different neutrino fla- 
vors. Thus the existence of transition magnetic moments 
implies that in presence of an external electromagnetic 
field conversions of left-handed neutrinos into their right- 
handed counterparts of another flavor (neutrino spin- 
flavor oscillations Q], 0) are possible, v~ -R> ft, where 
q/(3 and the index ± corresponds to different neutrino 
helicities. 

In this work we will consider spin-flavor oscillations of 
Dirac neutrinos in dense matter under the influence of 
strong external magnetic fields. Although the existence 
of Majorana neutrinos is favored in various neutrino mass 
generation models, like see-saw mechanism Q, the ques- 
tion whether neutrinos are Dirac or Majorana particles 
is still open 

Note that Dirac and Majorana neutrinos have com- 
pletely different structure of the magnetic moments ma- 
trix. Dirac neutrinos can have both diagonal and transi- 
tion magnetic moments whereas Majorana neutrinos only 
transition ones (see, e.g., Ref. @). Moreover the mag- 
netic moments matrix is symmetric in the Dirac case and 
anti-symmetric for Majorana neutrinos j5j. 

The strongest laboratory constraint on the neutrino ef- 
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fective magnetic moment, obtained by the GEMMA col- 
laboration [6j, is 3.2 x 10 -11 /ib, where hb is the Bohr 
magneton. Slightly weaker upper bound on the magnetic 
moments of neutrinos was obtained by the BOREXINO 
collaboration [7|. Note that the astrophysical limits on 
the magnetic moments of Dirac and Majorana neutri- 
nos extracted from the white dwarfs cooling rate Q can 
be even stronger than the laboratory ones. Therefore in 
most realistic situations, where the magnetic field is not 
too strong, neutrino spin and spin-flavor oscillations are 
likely to play a sub- leading role @ . Nevertheless in some 
astrophysical media the dynamics of neutrino oscillations 
can be significantly affected by a magnetic field because 
of its interaction with neutrino magnetic moments. For 
example, during a supernova explosion or in the vicinity 
of a neutron star, magnetic fields can reach values up 
to 10 16 G which is enough to influence the neutrino 
oscillations process in a significant way. 

Neutrino spin and spin-flavor oscillations in a super- 
nova were discussed in Refs. 

HEH. Th e impact of neu- 
trino magnetic moments and spin-flavor oscillations on 
r-process nucleosynthesis during a supernova explosion 
was studied in Ref. [12]. The neutrino spin flip, i.e. the 
transformation like v~ o 1/+ within the same flavor, can 
happen not only because of the magnetic moment inter- 
action with an external magnetic field, but also in colli- 
sions with background matter during a supernova explo- 
sion, that was examined in Ref. 13]. 

In this work we study neutrino spin-flavor oscillations 
in frames of the relativistic quantum mechanics [l4Ml8l j . 
Within the developed formalism we find the wave func- 
tions of neutrino mass eigenstates, which are the super- 
position of flavor neutrinos, for a given initial condition 
(Sec. [TTJ) . Using this method for the description of the 
neutrino evolution one can exactly take into account neu- 
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trino masses and the influence of external fields since 
we use exact solutions of the Dirac equation for mas- 
sive Dirac neutrinos. Note that neutrino oscillations in 
dense matter and strong magnetic fields can be also de- 
scribed with help of the methods of finite temperature 
field theory 

Then, in Sec. lIIIl we analyze the conventional quantum 
mechanical approach to the description of neutrino spin- 
flavor oscillations 0. We discuss the cases of the neutrino 
propagation in dense matter and weak magnetic field and 
in low density matter and strong magnetic field. Then 
we demonstrate that the relativistic quantum mechanics 
approach is consistent with the previously used formal- 
ism and obtain the corrections to the standard effective 
Hamiltonian. We show that the correction obtained re- 
sults in the appearance of the new resonance in neutrino 
oscillations. 

In Sec. IIVI we discuss a possible application of our re- 
sults to the conversion of right-handed electron neutrinos, 
which can be created in a supernova explosion [l3| , into a 
quasi-degenerate in mass sterile neutrino when particles 
propagate in an expanding envelope of a supernova and 
interact with its magnetic field. The importance of other 
factors, which can also influence the neutrino oscillations 
process, is considered in Sec. El Finally, in Sec. IVI[ we 
summarize our results. 



II. RELATIVISTIC QUANTUM MECHANICS 
DESCRIPTION 

Let us study the time evolution of the system of two 
mixed flavor neutrinos (y a ,vp) in matter and in an ex- 
ternal electromagnetic field F^ LU = (E, B). The indexes 
a and f3 can stand for any of the neutrino flavors e, jti or 
r. The interaction with background matter can be rep- 
resented in terms of the external axial vector fields f xx , . 
The Lagrangian for our system has the form, 

\=a,/3 \\'=a,P 

+ \M xx ,a flv F^y x ,, (2.1) 

where (toaa') is the nondiagonal mass matrix and {M\y) 
is the matrix of the neutrino magnetic moments. Note 
that in general case the matrices (m\\t) and (M\\i) are 
independent, i.e. the diagonal form of the matrix (m\\i) 
in a certain basis does not imply that of the matrix 
(M AV ). 

In the case of the standard model neutrino interac- 
tions the matrix (f xx> ) is diagonal: f xx , = f x 5\\>, where 
fx = (f x ,f\). The possible nondiagonal elements of the 
matrix f xx ,, with A ^ A', correspond to the non-standard 
neutrino interactions with matter [l7j . In the following 
we will study neutrino oscillations in non-moving and un- 
polarized matter with f\ = 0. The explicit form of the 
zero-th component /? = fx of the four vector f x for the 



electoneutral medium consistin g of electrons, protons and 
neutrons can be found in Ref. [2(| ■ In the case of sterile 
neutrinos f Us =0. 

To study the time evolution of the system (I2.1[) it is 
necessary to formulate the initial condition for the flavor 
neutrinos v\. We suppose that only one neutrino fla- 
vor, e.g., "j3", is present initially (see Refs. [T^ - [l8| ). i.e. 

u a (r, 0) = 0, and vp(r, 0) = 2/g (r), where (r) is the 
given function. Then we will look for the wave function 
v a at t > 0. If we choose a — fi or r and j3 = e, it will 
correspond to the typical situation of neutrinos emitted 
in the Sun: for the given initial flux of electron neutrinos 
we study the presence of other neutrino flavors at subse- 
quent moments of time. Note that the analogous initial 
condition problem for neutrino wave packets was studied 
in Ref. HJ. 

To analytically study the Dirac equation in presence of 
external fields we will consider the case of the coordinate 
independent functions f\ and F^. The analysis of the 
validity of this approximation will be made in Sec. [V] 
For the coordinate independent external fields the mo- 
mentum of the particles p is conserved. 

Moreover, the additional constraint can be imposed on 
the wave function i/^\r), 

P ± ^)(r)=^(r), ^=(l±^), (2-2) 

where £ = — j°j 5 ~f are the Dirac matrices. Eq. (|2.2j) 
implies that initially neutrinos of the flavor "/?" have a 
certain helicity. If we act with the operator on the 
final state v a (r, t) , we can study the appearance of the 
opposite helicity eigenstates among neutrinos of the fla- 
vor "a" , i.e this situation corresponds to the neutrino 
spin flavor oscillations 

Note that, unlike the helicity states defined in 
Eq. (|2.2|) . the chirality of a particle is a eigenvalue of 
the matrix 7 s : 7 5 ?/> = xij}. In case of massive Dirac par- 
ticles the helicity and the chirality do not coincide |22j . 
The helicity is a conserving number whereas the chiral- 
ity - not, since the operator 7 s does not commute with 
a Hamiltonian. 

Let us introduce the neutrino mass eigenstates ip a , a = 
1, 2, as 

a=l,2 v ' 

to diagonalize the mass matrix (rn\y). In Eq. (|2.3[) we 
take into account that for the two neutrinos system the 
mixing matrix (U\ a ) can be parameterized with help of 
one vacuum mixing angle 9. We suppose that the mass 
eigenstates tp a are Dirac particles with the masses m a . 

The Lagrangian (12.11) expressed in terms of the fields 
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ip a takes the form, 

a=l,2 

E ^{ll^ + ^uW^FA^, (2.4) 



26=1.2 



where 



(, ab )=uHM xx ,)u=^ £ 

{9a b )=UHh»)U=( 9l g ^ 



(2.5) 



are the nondiagonal matrices of neutrino magnetic mo- 
ments and neutrino interaction with matter in the mass 
eigenstates basis. In Eq. (|2.5[) we take into account that 
background matter is non-moving and unpolarized, i.e. 
only the zero-th component g a b = g^ b of the four vector 
9ab survives. 

To discuss the time evolution of the system (j2.4[) 
we write down the wave equations which result from 
Eq. dlU, 

iipa =H a ipa + Vip b , a = 1,2, a ^ b, 
U a =(ap) + f3m a - f i a f3J: 3 B + g a (l - 7 5 )/2, 
V = -nPZ 3 B + 3 (l- 7 5 )/2, (2.6) 

where a — 7°7 and /3 = 7 are the Dirac matri- 
ces. Here we study the neutrino motion along the x- 



axis: p = (p, 0,0), in only transversal magnetic field: 
B = (0, 0, B) and E = 0. 

Note that we cannot directly solve the wave equa- 
tions (12.6[) because of the nondiagonal interaction V 
which mixes different mass eigenstates. In vacuum, i.e 
in the absence of external fields, when (g^ b ) — and 
Ffj,u = 0, the mass eigenstates ^1,2 decouple and the sys- 
tem (|2.6I) can be easily solved. Nevertheless we can point 
out an exact solution of the wave equation iip a = Haipa, 
for a single mass eigenstate ip at that exactly accounts for 
the influence of the external fields. 



We look for the solution of Eq. 
form EMU: 



in the following 



Va{r,t)-e J ( 27r )3/2 e 



]T [aP(t)uPexp(-iEPt) 



C=±i 



(2.7) 



where the energy levels, which were found in Ref. (18j . 
have the form, 



(2.8) 



where R 2 a = yJp 2 M 2 + (p a B) 2 ml and M a = 
\J (/j, a B) 2 + g%jA. The basis spinors can be found in the 
limit of a small neutrino mass llal. 



,(0 



2y/2M a {M a - (ga/2) 



,(0 =. 



2^2M a (M a + Cg a /2) 



p a B + (M a - g a /2^ 
H a B - (M a + g a /2 
p a B - (M a + g a /2 

f l a B + (M a - ga/2 

( Ma-([HaB-g a /2] 

M a + C\f*aB + ga/2] 

-M a - C\MaB + g a /2] 

\-M a + CifiaB - g a /2] 



(2.9) 



It should be noted that the discrete quantum number £ = ±1 in Eqs. (|2.7|) - (|2.9|) does not correspond to the helicity 
quantum states. 

Now our goal is to find the time dependent coefficients aP(t) and bP(t). In the case of neutrino propagation in 
vacuum these coefficients do not depend on time and their values are completely defined by the initial condition. If we 
put the ansatz (|2.7j) in the wave equations (|2.6p . we get the following ordinary differential equations for the function 



-9b)t/2 



exp 



C'=±i 



.(C) 



exp 



-iEPt 



exp 



(2.10) 



To obtain Eq. (|2.10j) we use the orthonormality of the basis spinors (12. 9[) . Note that the differential equation for the 



function b a is analogous to Eq. (|2.10|) and thus omitted. Moreover, taking into account the fact that (u. 

we get that the eq 
energy eigenstates. 



(0 



V\v 



(C)\ 



= 0, 



we get that the equations for a a (t) and b a (t) decouple, i.e. the interaction V does not mix positive and negative 
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Let us rewrite Eq. (|2.10l) in the more conventional effective Hamiltonian form. For this purpose we introduce the 



"wave function" T = (a 1 , a 



2 J a i > a 2 



Directly from Eq. (|2.10j) we derive the equation for 



d*' 

i— = H'V, H' = 
dt 



( 


h- 


e iLU - t 














H+e -in+t 








H, 


e iSl+t 
















/ 



(2.11) 



where h± = (u£\V\u£), H± = (uf\V\uf), cu T = Ef - E* + ( 9l - g 2 )/2, and il T = Ef - Ef + ( 9l - g 2 )/2. Note 
that we do not give here the explicit form of the scalar products h± and H± in order not to encumber the text. 
Instead of ty' it is more convenient to use the transformed "wave function" defined by = U^f, where U = 

f2+)/2. Using the property lj + + w_ = 



diag{e i ( n +"-) t / 2 ,e i ( a -"-) t / 2 ,e- i (°-"+) t / 2 ,e- i ( n+ "+) t / 2 } and Q = 



+ -r we arrive to the new Schrodinger equation for the "wave function" 



d* 
l d^ 







h- 





H- 


-w_)/2 


H+ 





H+ 


-(il-u)+)/2 


h+ 





h+ 


-(n + oj + )/2 



(2.12) 



Despite initially we used the analog of the perturbation theory to analyze the influence of the potential V on the 
dynamics of the system (I2.6[) , the contribution of this potential is exactly taken into account in Eq. (|2.12l) . 

As we mentioned above, the quantum number £ does not correspond to a definite helicity eigenstate. Thus the 
initial condition, which we should add to Eq. (|2.12[) depends on the neutrino oscillations channel. Suppose that 



one has found the solution of Eq. (|2.12l) as \f rT (i) 
oscillations channel can be found as 



(V'li i } 2i "03, tpi)- Then the transition probability for Va 



•(*) 




VMT+gi/2 



y/Mi-gi/2 



y/M 2 + 92/2 y/M 2 - g 2 /2 



(2.13) 



To obtain Eq. (j2 . 13[) for simplicity we suppose that initially we have rather broad (in space) wave packet, corresponding 
to the initial condition ^''(r) ~ e lpr . 

Eqs. (|2.12p and (|2.13l) look completely new. However in Sec. lIIII we will show that the relativistic quantum mechanics 
approach to the description of neutrino spin flavor oscillations is consistent with the conventional quantum mechanical 
approach to this problem. 

I 

proach, developed in Sec. [TTl to the standard quantum 
mechanical treatment. 

Historically the Schrodinger equation which describes 
the dynamics of the Dirac neutrinos system was formu- 
lated in the flavor eigenstates basis 0, [4] . Definitely it is 
more convenient to use the flavor eigenstates basis since 
one gets the transition probability directly from the solu- 
tion of the Schodinger equation without additional ma- 
trix transformation (|2.3I) . Nevertheless we will formulate 
the dynamics of the neutrinos system in the mass eigen- 
states basis since the energies are well defined only for 
the neutrino mass eigenstates and one can distinguish 
the nature of neutrinos, i.e. say whether neutrinos are 
Dirac or Majorana particles, only in this basis. 

Let us discuss the situation when a neutrino moves 
in sufficiently dense matter and interacts with a weak 
magnetic field. The opposite case of a strong mag- 
netic field and a low density medium will be considered 
later. If the matter density is great, then in Eq. 



III. QUANTUM MECHANICAL DESCRIPTION 



In this section we analyze spin-flavor oscillations of 
Dirac neutrinos in frames of the standard quantum me- 
chanical approach. The main concept of this approach is 
the construction of an effective Hamiltonian acting in the 
space of quantum mechanical neutrino "wave functions" . 
Thus the proper choice of the basis of wave functions is 
as important as the correct form of the effective Hamil- 
tonian. In the majority of works devoted to neutrino 
spin-flavor oscillations the basis of helicity eigenstates 
was adopted (see, e.g., Ref. 0). As we will see below, 
this choice is justified only in the case of relatively weak 
external magnetic field and dense matter. We also con- 
sider the quantum mechanical derivation of the effective 
Hamiltonian in the opposite situation of strong magnetic 
field and low density matter. Finally we demonstrate 
the consistency of the relativistic quantum mechanics ap- 
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the averaged interaction of the diagonal magnetic mo- 
ment with an external magnetic field is less than the 
averaged diagonal interaction with background matter: 
(HaflYizB) -C (<7a(l — 7 5 )/2)- In this approximation it is 
convenient to rewrite the diagonal part of the Hamilto- 
nian in Eq. ()2.6[) as 



■Ha -^H a + V a , V Q = -//„/3E 3 B, 

H a =(ap)+/?m a +.g Q (l- 7 5 )/2, (3.1) 

to extract the small interaction V a from the main part 
of the diagonal Hamiltonian H a . We also rewrite the 
nondiagonal interaction V in Eq. (|2.6|) in the following 
form: V = Vb + V m , where Vb = —n(3T, 3 B and V m — 
g(l — 7 )/2, to separate the nondiagonal magnetic and 
matter interactions. 

One can notice that now the helicity operator de- 
fined in Eq. (|2.2p commutes with the modified Hamilto- 
nian H a - Therefore one can choose the helicity eigen- 
states basis for the quantum mechanical "wave func- 
tions" instead of the more complete set of basis func- 
tions (|2.9|) . If we study neutrinos propagating along 
the x-axis, these basis functions should satisfy the con- 
dition, (1/2)(1 ± £i)it± = ±u a z . The explicit form of 
these spinors can be found in Refs. [13, [Hj]. Note that 
we do not take into account the negative energy spinors 
Va since in Sec. [IT] we demonstrated that the evolution 
equations (|2.10[) for a a and b a decouple in matter and 
magnetic field. 

Now we construct the effective 4x4 Hamiltonian 
Hqm which governs the dynamics of the "wave func- 
tion" ^qIJ = (ipi , , ipf , 1P2 ) ) where ipf are the time 
dependent c-number wave functions corresponding to a 

definite helicity. The diagonal elements of Hq M can be 
calculated as the mean values of H a over the states with 
definite helicity Ua , i.e. they are equal to the energies 
of a neutrino moving only in background matter [l7l |23| , 



2p 



E + =p+^ + ..., (3.2) 



where we use the limit of ultrarelativistic particles. The 
nondiagonal elements of the effective Hamiltonian Hq M 
are the mean values of the operators V a , Vb and V m over 
the same helicity eigenstates Ua ■ Finally we arrive to 
the effective Hamiltonian in the ultrarelativistic limit, 



n QM — 



$ + 91 


9 




-fj,B~ 


9 


-$ + 92 


-fiB 


-IMtB 


-H X B 


-fxB 


$ 















(3.3) 



where $ = dm 2 /4k is the phase of vacuum oscillations 
and dm 2 — m? — m 2 ,. The Hamiltonian H^\, determines 



QM 



(i) 

QM- 



the time evolution of the "wave function" ^> 

One can see that we have reproduced the standard 
effective Hamiltonian proposed in Ref. to study spin- 
flavor oscillations of Dirac neutrinos. In our analysis it 



is important that the helicity eigenstates u a are used 
as the basis functions. It is correct only if the diagonal 
magnetic interaction (fi a f3T,^B) is small. The numerous 
works devoted to neutrino spin- flavor oscillations (see, 
e.g., the recent review 24]) showed that the Hamilto- 
nian (|3.3[) seems to be applicable to a more general situ- 
ation. Nevertheless here we show that the derivation of 
this Hamiltonian presented, e.g., in Ref. Q is restricted 
to the case of small diagonal magnetic interaction. 

Now we discuss the opposite case when a neutrino in- 
teracts with a very strong magnetic field and moves in low 
density matter. Hence in Eq. (|2.6[) the diagonal magnetic 
interaction is much bigger than the diagonal interaction 
with background matter: (fx a ^T,sB) ^> (g a (l — 7 5 )/2). 
In this case it is also convenient to redefine the diagonal 
Hamiltonian in Eq. (|2.6p in the following way: 



Ha^Ha+Va, V a - , 9a (1 - 7 5 )/2, 

H a =(ap) + /3m a - fi a ^ 3 B. 



(3.4) 



Note that now the Hamiltonian H a does not conserve the 
helicity of a particle. Thus to develop the standard quan- 
tum mechanical approach we have to choose the proper 
basis in the space of neutrino wave functions. For the ba- 
sis wave functions it is convenient to use the eigenvectors 
of the operator [16|, n a = m a S3 + i7°7 5 (S x p) 3 — [i a B, 
which characterizes the spin direction with respect to the 
magnetic field. The explicit form of these spinors is given 
in Ref. 0. 

The effective 4x4 Hamiltonian Hq^ acting in the 
space of the quantum mechanical neutrino "wave func- 
tions" = (i'l 1^2 i^Pi ^t) can be constructed in 
a straightforward way as in the case of the weak diag- 
onal magnetic interaction. Here tp^ are the c-numbcr 
time dependent functions representing neutrino states 
with a definite spin projection on the magnetic field 
direction. However we can notice that the Hamilto- 

(2) 

nian Hq M can be obtained by the similarity transfor- 
mation, HqIj = V2iHqI j VJ 1 , with the orthogonal ma- 
trix V21 = (7 5 — 7°)/v / 2, with the Dirac matrices 7 and 
7 5 taken in the standard representation [25]. It should 

be noted that the direct calculation of the Hamiltonian 

(2) 

HqIj gives the same result. 

Despite we demonstrated the consistency of the de- 
scriptions of neutrino oscillations in different bases there 
is a conceptual discrepancy between these cases. The 
quantum mechanical treatment of neutrino evolution re- 
quires the correct choice of both the effective Hamiltonian 
and the basis of "wave functions" where this Hamiltonian 

(1 2)T 

acts. The quantum mechanical "wave functions" ^q'm 
are applicable to the description of neutrino evolution 
for the cases of weak and strong diagonal magnetic in- 
teraction respectively since they correspond to different 
conserved quantum numbers in each case. Although for 
ultrarelativistic neutrinos the description of neutrino os- 
cillations is identical in these bases, there can be a differ- 
ence if we calculate the corrections to the leading term. 
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Now let us check the consistency between the relativistic quantum mechanics method developed in Sec. Ql] and the 
standard quantum mechanical description of neutrino oscillations. We will study ultrarelativistic neutrinos. In this 
approximation the energy levels (|2.8j) take the form, 



2 2 2 

2 2p 8p 2 M a 



(3.5) 



where we keep the term ~ m 2 /k 2 to examine the corrections to the leading term. 

Performing the similarity transformation of the effective Hamiltonian H in Eq. (|2.12[) with the orthogonal matrix 
V of the following form: 



V 



f-^/l + gi/2M 1 

J_ -^/l + .92/2X2 

V2 V 1 -5i/2Mi 

V y/l - g 2 /2M 2 



\ 



l i 1 B/^M 1 (M 1 +g 1 /2) __0 

^2S/ v /X 2 (X 2 + . 92 /2) 

HiB/y/Mi(Mi-gi/2) 

^B/^M 2 {M2-92/2) J 



(3.6) 



we can see that the Hamiltonian transforms to V T HV ~ Hqm + 8H, where Hqm — Hqm 
the 4x4 unit matrix, and 



tr{H$l)/4 ■ I, I is 



SH 



1 



16k 2 



diag 



2 9i 



2 fl 2 



2 5i 



2 52 



" 2 M 2 2 ' mi Mt' m2 Ml 



(3.7) 



is the correction to the standard effective Hamiltonian. 
It should be noted that the transformation matrix V in 
Eq. p.6j) depends on the magnetic field strength and the 
matter density, whereas the matrix V21 is external field 
independent. 

The effective Hamiltonian Hqm is equivalent to Hgj^ 
in Eq. Q3.3P since the unit matrix does not change the 
particles dynamics. If we omit the correction p.7j) . that 
is valid for ultrarelativistic neutrinos, we can see that the 
relativistic quantum mechanics approach is equivalent to 
the standard quantum mechanical method. The correc- 
tion 8H results from the fact that in Sec. [II] we use the 
correct energy levels (|2.8p and the wave functions (|2.9p 
for a neutrino moving in dense matter and strong mag- 
netic field. 

Note that in Eq. (|3.7p we keep only the diagonal correc- 
tions to the effective Hamiltonian (|3 .3[) since the appear- 
ance of a resonance in neutrino oscillations is sensitive 
to the diagonal elements a Hamiltonian. We should re- 
mind that the expressions for the basis spinors (|2.9p were 
obtained in the approximation of small masses of neutri- 
nos, whereas in the expansion of the energy levels Q3.5P 
we keep terms up to ~ m 2 /k 2 . However using Eqs. (|2.12|) 
and (|3.6I) we get that the additional small contributions 
~ m 2 /k 2 from the basis spinors are washed out in diag- 
onal entries in Eq. (|3.7p . 

IV. APPLICATIONS 

In this section we study the application of the ob- 
tained Hamiltonians (|3.3p and (13.71) to the description 
of oscillations between electron and hypothetical ster- 
ile neutrinos in an expanding envelope of a supernova. 



The existence of sterile neutrinos closely degenerate in 
mass with electron, muon- or tau-neutrinos was recently 
discussed in Ref. [26| in connection to solar and super- 
nova neutrinos. The mass squared differences consid- 
ered in these publications were in the following range: 
I0" 19 eV 2 < 5m 2 < I(T 12 cV 2 . 

On the contrary, the mixing angles of these additional 
neutrinos cannot be well constraint. Therefore can as- 
sume that vacuum mixing angle is small, 9 -C 1. In 
this case mixing matrix between mass ip a) a = 1, . . . , 4, 
and flavor v\, A = e, /i, r, s, neutrino eigenstates has the 
form U4 sa diag ([/3, 1), where U3 is the mixing matrix 
of the three Dirac neutrinos system (see Eq. (|2.3p and 
Ref. [27|). It should be noted that it is very difficult 
to detect additional neutrino flavor if it is weakly mixed 
with active neutrinos and Sm 2 is small. Such a neutrino 
can be revealed through spin-flavor oscillations only if it 
has a transition magnetic moment. 

Besides the huge amount of left-handed neutrinos from 
a supernova, a smaller flux of right-handed particles is 
predicted [I3( . These right-handed neutrinos can be cre- 
ated in the following reaction: + (e~,p,N) — >• v + + 
(e~ ,p, N), with electrons e~, protons p, and nuclei N in 
the dense matter of a protoneutron star. For the neu- 
trino spin-flip in matter to happen within one generation 
v~ — > a neutrino should be the Dirac particle with 
a nonzero diagonal magnetic moment. Left-polarized 
supernova neutrinos are strongly degenerate and can 
occupy energy levels above the Fermi surface. Right- 
polarized neutrinos created in the spin-flip reactions have 
the energy of the order of the left-polarized neutrinos. 
The detailed analysis of Ref. shows that the energy 
of these particles is in the E v — (100 — 200) MeV range. 

The generation of electron neutrinos with right-handed 
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polarization can provide additional supernova cooling 
since they do not interact with background matter and 
thus freely carry away the supernova energy. Moreover 
these particles can be potentially detected in a terres- 
trial detector due to their spin precession, back to left- 
handed states, in the galactic magnetic field. However, if 
we point out an additional neutrino oscillations channel, 
which contributes to the right-handed neutrinos dynam- 
ics, these particles are unlikely to be observed. 

Let us study the appearance of a resonance in v+ «-» v~ 
oscillations channel. It is known that a resonance in a 
certain channel of neutrino oscillations can appear if the 
difference between two diagonal elements in the effective 
Hamiltonian is small [28[. Using Eqs. (|3.3|) and Q3.7P 
in the approximation of 8 <C 1, as well as the results 
of Ref. [2(| we obtain the resonance condition for these 
oscillations as, 

Sm 2 «5.0 x 10~ 17 eV 2 x (3Y e - 1) 




where Y e — n e /(n e + n n ) is the electrons fraction, p is 
the mass density of background matter, and m„ e is the 
absolute mass of an electron neutrino. Since 5m 2 is sup- 
posed to be small, we can set m v to be equal to either 
mi or the mass the additional neutrino mass eigenstate 
appearing due to the presence of a sterile neutrino. In 
Eq. (|4.1I) we suppose that matter is electroneutral and 
the diagonal magnetic moment of an electron neutrino is 
small: p Ve B <C / e . 

Suppose that the flux of right-handed electron neutri- 
nos is crossing an expanding envelope of a supernova. A 
shock wave is usually formed in the envelope 29]. Ap- 
proximately 1 s after the core collapse, the matter den- 
sity in the shock wave region L ~ 10 8 cm can be up to 
10 6 g/cm 3 . We can also suppose that the matter den- 
sity is approximately constant inside the shock wave. 
From Eq. (|4.1[) we can see that a resonance in neu- 
trino oscillations happens if p ~ 10 6 g/cm 3 , Y e > 1/3, 
E v ~ 100 MeV, m Ue « 0.2 eV, and <5to 2 « 2 x 10~ 18 eV 2 , 
that is consistent with the modern cosmological limits 
on the absolute neutrino mass [3(| , the mentioned above 
estimates of the mass squared differences with a sterile 
neutrino and the energy of right-handed neutrinos, as 
well as to the parameters of a shock wave. 

Besides the fulfillment of the resonance condition (|4.1I) , 
to have the significant v+ O v~ transitions rate the os- 
cillations length should be comparable to the the shock 
wave size L. We can express this condition in the follow- 
ing form: 

B«5.3x 10 7 Gx f ^ ^ (— £— ) -(4-2) 

\IQ- 12 p B J V^kmy v ' 

For the transition magnetic moment p — 3 x 10 -12 pb [HI 
and L ~ 10 8 cm (see above), we get that B ~ 10 7 G. 
Supposing that the magnetic field of a protoneutron star 



depends on the distance as £?di p (r) = Bo(R/r) 3 , where 
R = 10 km is the typical protoneutron star radius and 
Bo = 10 13 G is the magnetic field on the surface of a pro- 
toneutron star, we get that at r = 10 s cm the magnetic 
field reaches 10 7 G, which is consistent with the estimates 
of Eq. (|4.2[) . Note that magnetic fields in a supernova ex- 
plosion can be even higher than 10 13 G and reaches the 
values of - 10 16 G 0. 

In Fig. []Ja) we present the numerical solution of 
the Schrodinger equation with the Hamiltonians (|3.3[) 
and (|3.7[) taking into account the diagonal magnetic 
moment of an electron neutrino which was neglected 
while obtaining the resonance conditions since we as- 
sumed that p a <C p. It is necessary to remind that in 
Ref. [l3| it was found that to get a significant luminos- 
ity — 10 50 erg/s the diagonal magnetic moment should 
be p Vc — 10~ 13 ^ B . We suppose that both Eqs. (|4.1j) 
and (|4.2[) are satisfied. One can see in Fig. HJa) that di- 
agonal magnetic moments do not significantly influence 
the dynamics of spin-flavor oscillations since the numer- 
ical transition probability practically coincide with the 
approximate analytical expression P(x) — sin 2 (pBx) ob- 
tained from Eqs. ([3~3]), ([377]), (gTT]), and g~2]) in the limit 
of small p Ue . 

We plot Fig. [IJa) in the approximation of constant 
external fields. Despite the matter density decreases as 
1/r 3 inside the envelope (32j, we can suppose that it is 
approximately constant in the shock wave region. How- 
ever the approximation of the constant magnetic field is 
quite rough. In Fig. |Tfb) we present the solution of the 
exact Schrodinger equation for the coordinate dependent 
magnetic field Bdi P (r), for various values of the magnetic 
fields Bq at the neutrinosphere. We again suppose that 
the resonance condition (I4.ip is fulfilled. 

Note that in Sec. |H] we assumed that the magnetic 
field is coordinate independent. However, if we suggest 
that the length of the spatial variation of a magnetic field 
is much bigger than the typical size of a neutrino wave 
packet (see also Sec. |Vj, we can neglect it in the anal- 
ysis of the neutrino wave equation (|2.6p . Thus, if this 
condition is satisfied, the magnetic field is supposed to 
be spatially constant in deriving main Eqs. (|2.12[) , Q3.3p . 
and ([377]) in Sees. M and Hill Nevertheless we take into 
account the magnetic field variation while solving the ef- 
fective Schrodinger equation with the Hamiltonians (|3.3p 
and (|3.7p . This requirement for spatial variations of mag- 
netic fields can be characterized as the microscopic adi- 
abaticity condition in contrast to the macroscopic adia- 
baticity typically used in the analysis of neutrino oscilla- 
tions (see, e.g., Ref. |33jV 

One can see in Fig. QJb) that at relatively weak mag- 
netic fields Bq ~ 10 13 G the transition probability is a 
strictly increasing function reaching the asymptotic value 
~ 0.5 and never is equal to one. This magnetic field corre- 
sponds to B(r — i?cnv) = 10 7 G, where R env = 10 8 cm is 
the internal radius of an expanding envelope. We should 
remind that oscillations in a constant magnetic field of 
such a strength are at resonance, cf. Eq. (|4.2I) . and hence 
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FIG. 1. (a) The numerical transition probability versus the distance traveled by a neutrino beam when both resonance 
conditions, Eqs. ()4.1[) and H4.2|l . are satisfied. The parameters of the neutrino system are the following: n = 3 X 10 _12 /Ub, 
\j, Ve — 10~ 13 /ib, fJ-v s = 0, E — 100 MeV, m Ve = leV and 8m 2 = 5 x 10 _17 eV 2 . Neutrinos move in background matter with 
p = 10 6 g/cm 3 and interact with the magnetic field B = 10 7 G. (b) The numerical transition probability as a function of the 
distance in the case of the coordinate dependent magnetic field -Bdip(f) for various values of Bo- The distance x is measured 
from the inner radius of an envelope R C nv = 10 s cm. The parameters of the neutrino system are the same as in the panel (a) . 



the transition probability can reach a unit value. 

We also notice that at the big distances traveled by 
neutrinos the transition probability becomes constant. 
Indeed, using Eqs. (I3.3[) and Q3.7P we get that in the 
limit of the small diagonal magnetic moment of an elec- 
tron neutrino and if the resonance condition (I4.ip is ful- 
filled, at big distances the transition probability has the 
form, P(x > R cnv ) = sin 2 (^B -R 3 /2i?onv)- Tne values of 
the asymptotic transition probability shown in Fig. [1] are 
in agreement with this expression. At strong magnetic 
fields (B = 10 14 G or B{r = R cnv ) = 10 s G) the tran- 
sition probability can reach big values at the outer edge 
of a broad envelope, r ~ 10 9 cm. Note that analogous 
behaviour of transition probability was found in Ref. [34fl 
while studying spin-flavor oscillations of Majorana neu- 
trinos in a supernova. 



V. ANALYSIS OF APPROXIMATIONS 

First we should remind that we use the relativistic 
quantum mechanics approach, with the external fields be- 
ing independent of spatial coordinates. If external fields 
depend on the spatial coordinates, Dirac wave packets 
theory reveals various additional phenomena such as par- 
ticles creation by the external field inhomogeneity [35| . 
For the approximation of the spatially constant exter- 
nal fields to be valid, the typical length scale of the ex- 
ternal field variation L cxt should be much greater than 
the Compton length of a neutrino [35|: L cxt >• Xc — 
h/m v c 20]. For a neutrino with m„ ~ leV this condi- 
tion reads L ext ^S> 10~ 5 cm, that is fulfilled for almost all 
realistic external fields. 

We should also make a remark on the accounting for 



the mixing potential V in Eq. (|2.10p . Eq. (|2.6p contains 
the nondiagonal term ~ V. The analysis of analogous 
equations is typically made in frames of the perturbation 
theory (3(| using the expansion on powers of the mixing 
potential, i.e. on the coupling constants which are pro- 
portional to g and p:B. As it was mentioned in Sec. [TTl 
the wave equations for massive neutrinos in vacuum de- 
couple and the evolution of these states depends on the 
initial condition only. While solving Eq. (12.10)) we could 
have taken into account the terms linear in g and /iB 
as it was made in Refs. [lH [Iff]. However in subsequent 
calculations, which lead to Eqs. (|2.1ip and (|2.12l) . these 
coefficients were accounted for exactly. 

Now let us discuss other factors which can also give 
the contributions, comparable with Eq. (13.71) . to the ef- 
fective Hamiltonian. While deriving the effective Hamil- 
tonian (|2 . 12[) in Sec. |TT] we supposed that the magnetic 
field is transverse with respect to the neutrino motion. 
The effect of the longitudinal magnetic field B\\ on neu- 
trino oscillations was studied in Ref. [37[. In order to 
neglect the longitudinal magnetic field contribution in 
comparison with our correction (13.71) , its strength should 
satisfy the condition, 



Bu 



B± lQkB±\p: a m a - fi b m b \ 



m l9l m b9b 



Ml 



Ml 



, (5.1) 



where B± is the transverse component of the magnetic 
field. The condition (|5.ip is satisfied for neutrinos emit- 
ted inside the solid angle near the equatorial plane with 
the spread 2d ~ 2B\\/B^. Assuming the radially sym- 
metric neutrino emission we find that about 2% of the to- 
tal neutrino flux is affected by the new resonance (14.11) . 
i.e. the influence of the longitudinal magnetic field is 
negligible for oscillations of such particles. 
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The next important approximation made in the de- 
viation of Eq. (|2.12[) was the assumption of negligible 
polarization of matter which can be not true if we study 
rather strong magnetic fields. The effect of matter polar- 
ization on neutrino oscillations was previously discussed 
in Refs. [2(], HH, [3!|. It was found that in the lead- 
ing order in m a /k matter polarization produces the fol- 
lowing contribution to the diagonal entry of the effec- 
tive Hamiltonian, corresponding to right-handed neutri- 
nos: g a {\f(3 v ){m a /k), where /3„ is the neutrino velocity 
and A f is the mean polarization vector of background 
fermions. Note that in polarized matter the effective en- 
ergy of left-handed neutrinos also changes. However this 
process does not contribute to the considered oscillations 
channel v+ -R* v~ . 

It is clear that one should take into account only the 
polarization of electrons since nucleons are much heavier. 
Using the results of Refs. [HI, H(| we obtain that in the 
shock wave region electrons are relativistic and weakly 
degenerate. Finally we get that the new correction to 
the effective Hamiltonian (|3 . T[) becomes bigger than the 
contribution of matter polarization to the effective poten- 
tial of right-handed neutrinos if the electron temperature 
exceeds 10 MeV which is comparable with the mean tem- 
perature in an expanding envelope |40j . 

The presence of relatively big Dirac neutrino magnetic 
moments implies the existence of new interactions, be- 
yond the standard model, which electromagnetically cou- 
ple left- and right-handed neutrinos. It is possible that 
these new interactions also contribute to the effective 
potential of the right-handed neutrino interaction with 
background matter. Despite this additional effective po- 
tential is likely to be small, one should evaluate it and 
compare with the correction (|3.7|) . 

The most generic SU(2) L x U(l)y gauge invariant and 
remormalizable interaction which produces Dirac neu- 
trino magnetic moment was discussed in Ref. 41]. The 
effective Lagrangian of this interaction involves the di- 
mension n = 6 operators Oj, C c s — (1/A 2 ) J2j CjOj + 
h.c, where A ~ 1 TeV is scale of the new physics, Cj are 
the effective operator coupling constants and sum spans 
all the operators of the given dimension. 

One of the operators Oj also contributes to the ef- 
fective potential of a right-handed neutrino in matter, 
O = nLT a 4>a^ v vi{W^ v , where n is the coupling con- 
stant, T a are Pauli matrices, L T = (y^, is the SU(2)^ 
isodoublet, <f> — ir 2 0*, with <p being a Higgs field, and 
= - d v W£ - Ke abc W*WZ is the SU(2) L field 

strength tensor. Assuming the spontaneous symmetry 
breaking at the electroweak scale, (f> T — > (0,v/^/2), we 
can rewrite the contribution of the operator O to the 
effective Lagrangian in the form, 

£ off = ^-e L a^v R (W^ - iW* v ) + h.c, (5.2) 

which implies that a process e~ + v + — > e~ + v_ should 
happen in background matter. 



Using the results of Ref. [4l| we can evaluate the con- 
tribution of the new interactions to the effective Hamil- 
tonian (|3.3j) as 

G F M^ \ /ib / \m e J 

where V sm ~ Gpn e is the standard model effective po- 
tential, Gf is the Fermi constant, and Mw is the W 
boson mass. Taking fi u ~ 10~ 12 /is, E v ~ 100 MeV and 
m„ ~ 0.2 eV (see Sec. IIV|) we can get that the ratio of the 
correction to the effective potential SVr and new correc- 
tion (|3.7p is ~ 10~ 2 . It means that the influence of new 
interactions, which generate neutrino magnetic moments, 
are not important for neutrino spin-flavor oscillations. 

The constraint on the Dirac neutrino magnetic mo- 
ment obtained in Ref. [4l[ is 10 _14 /zb- Neverthe- 
less in Sec. IIVI we used the magnetic moments in the 
(10~ 13 — 10~ 12 )hb range since analogous constrains on 
the Dirac neutrino magnetic moments were obtained in 
Refs. [HI, [3l| on the basis of the analysis of astrophysical 
data. 



VI. CONCLUSION 

In this paper we have described neutrino spin-flavor 
oscillations in dense matter and strong magnetic field in 
the frame of relativistic quantum mechanics. The advan- 
tage of this formalism, compared to the commonly used 
quantum mechanical approach, is that one can exactly 
take into account the neutrino properties like initial mo- 
mentum, dm 2 , mixing angles and magnetic moments, as 
well as the matter density and the strength of the mag- 
netic field since we used the exact solutions of the Dirac 
equation for massive neutrinos in presence of external 
fields. 

In Sec. |TT] it was demonstrated that the initial condi- 
tion problem for the system of two mixed flavor neutri- 
nos, each of them represented as four-component Dirac 
spinors, can be reduced to a Schrodinger like equation 
[see Eq. (|2.12|) ]. In Sec. Illll it was shown that the Hamil- 
tonian of this evolution equation formally coincides with 
the previously proposed [2| effective Hamiltonian (|3 . 3[) . 

It should be, however, noted that the dynamics of neu- 
trino spin-flavor oscillation in matter and magnetic field 
in frames of the quantum mechanical description is de- 
fined by both the effective Hamiltonian and the correct 
basis of the neutrinos wave functions. In Sec. IIIII we 
showed that the choice of the helicity eigenstates basis is 
justified only in the case of the weak magnetic field limit. 
This kind of basis was used in the previous treatment of 
neutrino spin-flavor oscillations Q. The opposite situa- 
tion of the strong magnetic field and low density matter 
was also discussed in Sec. IIIII In this limit the correct 
basis consists of the eigenfunctions of the operator n a . 

Besides the demonstration of the consistency of our 
results with the standard approach for the description 
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of neutrino spin-flavor oscillations, we found the cor- 
rection (13.71) to the commonly used effective Hamilto- 
nian (|3.3|) which is usually omitted 0. It was possible 
to obtain this correction in the explicit form since we 
used the energy levels (12 .8[) and basis spinors (12 .9[) which 
exactly account for external fields. 

In Sec. lIVI we discussed the realistic situation when the 
correction Q3.7P is important. We considered i/+ <-» v~ 
oscillations channel, where v s is the additional ster- 
ile neutrino almost degenerate in mass with other neu- 
trino states. The right-handed electron neutrinos were 
supposed to be produced during the supernova explo- 
sion due to the scattering of left-handed neutrinos with 
the non-zero diagonal magnetic moment on background 
fermions [l3j . The flux of these right-handed neutrinos 
was taken to propagate through the expanding envelope 
and interact with an external magnetic field due to the 
presence of the non-zero transition magnetic moment. 

We found that the new resonance in neutrino spin- 
flavor oscillations can appear if the strength of the mag- 
netic field and the matter density are close to the values 
recently discussed in Refs. [l(| Note that new res- 
onance condition (|4.1I) depends on the absolute value of 
the neutrino mass. Thus the observation of this new res- 
onance would provide the information about the absolute 
scale of the neutrino masses. Although i/+ f-> v~ oscilla- 
tions do not change the dynamics of a supernova explo- 
sion the appearance of this additional resonance chan- 
nel of neutrino spin-flavor oscillations makes impossible 
a terrestrial observation of the right-handed supernova 
neutrinos proposed in Ref. [HJ]. 

To analyze the dynamics of the Schrodinger equation 
with exact Hamiltonians (|3.3p and (|3.7|) in Sec. IIVI we 
presented the numerical transition probability which ac- 



counts for all magnetic moments and is built for vari- 
ous magnetic fields configurations. In particular we an- 
alyze the constant magnetic field and more realistic co- 
ordinate dependent magnetic field of a magnetic dipole 
-BdipW- ^ was revealed that diagonal magnetic moment 
of an electron neutrino does not significantly influence 
the transition probability. Then we found that at rather 
strong magnetic field strength the transition probability 
can reach big values when a neutrino leaves an expanding 
envelope. 

In Sec. [V] we considered other factors which can be 
comparable with the correction obtained (|3.7[) . We an- 
alyzed the contributions of the longitudinal magnetic 
field and matter polarization which were omitted in the 
derivation of the effective Hamiltonian (|2.12[) . In partic- 
ular we have found that a longitudinal magnetic field is 
not important for neutrinos emitted near the equator of 
a star. The matter polarization does not influence the 
dynamics of neutrino oscillations if the temperature of 
background electrons is higher than a few MeV. We also 
evaluated the possible contributions of the new interac- 
tions, which generate neutrino magnetic moments [4lj . to 
the effective potentials of right-handed neutrinos. It was 
found that these contributions are negligible compared 
to the correction fl3.7[) . 
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